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Abstract 

The problem of cosmological particle creation for a spatially flat, homogeneous and isotropic 
Universes is discussed in the context of f(R) theories of gravity. Different from cosmological 
models based on general relativity theory, it is found that a conformal invariant metric does not 
forbid the creation of massless particles during the early stages (radiation era) of the Universe. 
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I. INTRODUCTION 



It is widely believed that the possible emergence of space and time trough a cosmic 
singularity and the presence of horizons in the Friedman-Robertson- Walker (FRW) models 
strongly suggest that matter and radiation need to be somehow created in order to overcome 
some basic conceptual difficulties of big-bang cosmology [l| . In these connections, the process 
of matter creation in an expanding universe has been extensively discussed in the last four 
decades either from a macroscopic, as well as from microscopic viewpoints [2-18]. 

Macroscopically, the matter creation was extensively investigated as a byproduct of bulk 
viscosity processes near the Planck era as well as in the slow-rollover phase of the new infla- 
tionary scenario , [5|. Later on, the first self-consistent macroscopic formulation of the 



matter creation process was put forward by Prigogine and coworkers 



6| and formulated in 



a covariant way by Calvao et al. 3| with basis on the relativistic nonequilibrium thermody- 
namics. In comparison to the standard equilibrium equations, the process of creation at the 
expense of the gravitational field is described by two new ingredients: a balance equation for 
the particle number density and a negative pressure term in the stress tensor. Such quanti- 
ties are related to each other in a very definite way by the second law of thermodynamics. 
In particular, the creation pressure depends on the creation rate and may operate, at level of 
Einstein's equations, to prevent either a spacetime singularity js Is] or to generate an early 



ain the 



inflationary phase [9|. More recently, such formulation has also been applied to exp 
late time accelerating stage of the Universe and other complementary observations [id, 111 ]. 

Microscopically, after the pioneering work of Parker 12] , the quantum process of particle 
creation in the course of the cosmological expansion has also been studied by several authors 
B Q, H, lis 113, ll8( . These Parker Particles' are created because, according to the co- 



variant equations of the fields in the Heisenberg picture, the positive and negative frequency 
parts of the fields become mixed during the universe expansion, so that the creation and 
annihilation operators at one time t\ are linear combinations of those ones at an earlier time 
t2, resulting at first, in a particle production. 

One of the most interesting results from Parker's work is that in a radiation dominated 
Universe, the positive- and negative- frequency mode functions for massless fields are not 
mixed in the course of the expansion due to the conformal invariance of the metric. In 
particular, this means that there is no creation of massless particles, either of zero or non- 
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zero spin. Viewed in another way, the dispersion relation in Fourier space is exactly the same 
one of the flat Minkowski space, and, as such, a wave propagation of a massless particle is not 
accompanied by particle production. As a consequence, in the GR framework, no photon, 
graviton or any other kind of massless particles are produced by purely expanding effects at 
early times [17| . 

Usually, the calculations of particle production deal with comparing the particle number 
at asymptotically early and late times, or with respect to the vacuum states defined in two 
different frames and do not involve any loop calculation. However, the main problem one 
encounters when treating quantization in expanding backgrounds concerns the interpretation 
of the field theory in terms of particles. The absence of Poincare group symmetry in curved 
space-time leads to the problem of the definition of particles and vacuum states. The problem 
may be solved by using the method of the diagonalization of instantaneous Hamiltonian by 
a Bogoliubov transformation, which leads to finite results for the number of created particles 

On the other hand, the major developments concerning the study of particle produc- 
tion in curved background are concentrated in the standard general relativity. Recently, 
non-standard gravity theories have been proposed as an alternative to explain the present 
accelerating stage of the universe only with cold dark matter, that is, with no appealing to 
the existence of dark energy. This reduction of the so-called dark sector is naturally obtained 
in f(R) gravity theories. In such an approach, the curvature scalar R in the Einstein- Hilbert 
action is replaced by a g eneral function f(R), so that the Einstein field equation is recovered 
as a particular case [20|, |21] . 

A basic motivation for this kind of generalization is the fact that higher order terms in 
curvature invariants have to be added to the effective Lagrangian of the gravitational field 



when quantum corrections are taken into account [22|, |23j, l24j . Two different formalisms 
are usually adopted to study f(R) theories. In the so-called metric approach, the equations 
of motion are obtained by varying the metric tensor, while in the Palatini formalism, one 
considers the metric and the affine connection to be independent of each other, and one 
has to vary the action with respect to both of them. These two methods lead to the same 
equations only if f(R) is a linear function of R (GR case). The metric approach leads to 
fourth-order equations driving the evolution of the scale factor while the Palatini method 
leads to second-order equations, so that it is appealing because of its simplicity. Actually, for 
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many interesting examples, it has been argued that some problems of instability are present 



in the metric approach 
such instabilities 



25j | , although some recent works have cast doubts on the existence of 



261 ] . In what follows, we will focus our attention on the Palatini formalism 
and a simple power law type of f(R)- As shown by Vollick 2^] using the Palatine approach, 
the inclusion of a term in the action leads to a theory of gravity that predicts late 
time accelerated expansion of the Universe. Sotiriou [27] also demonstrated that a R 3 term 
can account for an early time inflation, contrary to a R 2 type discussed by many authors 
(see Meng and Wang [28j and Refs. there in). More recently, Miranda et al. discussed a 
viable f{R) cosmology based on a logarithm contribution of the scale factor which emerges 
naturally as a limiting case of a general power law 29 1. 

In this letter we will present only the basic ideas and results on the scalar particle creation 
on expanding universe in the view of a R n type f{R) theory. In brief, our basic result can 
be stated as follows: massless particles in a radiation dominated Universe can be produced 
by purely expanding effects in generic f{R) theories. In other words, Parker's result is valid 
only in the context of general relativity. 



II. QUANTIZATION IN EXPANDING BACKGROUNDS 

The canonical quantization of a scalar field in curved backgrounds follows in straight 
analogy with the quantization in a flat Minkowski background with the gravitational metric 
treated as a classical external field which is generally nonhomogeneous and nonstationary. 
Let us summarize the basic results for a linear Klein-Gordon scalar field in a spatially flat 
FRW geometry (in our units c = 1). 

A real minimally coupled massive scalar field <fr is described in a curved spacetime by the 
action 



S = - I v ' —gd^x 



g a0 d a <pd <p - m 2 2 



2 , v . - (1) 

In terms of the conformal time T], the metric tensor g^ v is conformally equivalent to the 
Minkowski metric rj^, so that the line element is ds 2 = a 2 (r))r]^ l/ dx tJ 'dx u , where a(rj) is the 
cosmological scale factor. Writing the field (j)(rj,x) = a(r])~ 1 x, the equation of motion that 
follows from the action ([[]) is 



a" 



x"-V 2 x+[m 2 a 2 --) X = 0, (2) 
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where the prime denotes derivatives with respect to the conformal time 77. We can see 
that the field \ obeys the same equation of motion as a massive scalar field in Minkowski 
space-time, but now with a time dependent effective mass, 



a" 



[ ff ( v )=m 2 a 2 --. (3) 



This time dependent mass accounts for the interaction between the scalar field and the 
gravitational field. The energy of the field x is n °t conserved (its action is explicitly time- 
dependent), and, more important, its quantization leads to particle creation at the expense 
of the classical gravitational background. However, for a FRW type Universe dominated by 
radiation one finds that a" = 0. Therefore, for a massless scalar field, there is no particle 
production since Eq. (|2j) reduces to same one of the Minkowski spacetime. This is the Parker 



result which was clearly deduced with basis on the general relativity [121 ] . In other words, 
by assuming the general relativity, the gravitational field of the Universe at early stages is 
unable to work as a pump supplying energy for massless scalar fields. 

After a Fourier expansion of Eq. (j2J), we are left with the mode equation 

Xk + "l(v)Xk = 0, (4) 

with ujI(t]) = k 2 + m 2 [ jj-. The solutions of the above equation give the positive- and negative- 
frequency modes. Following standard lines, the quantization can be carried out by imposing 
equal-time commutation relations for the scalar field \ and its canonically conjugate mo- 
mentum 71 = x', namely [x{%, v) > ^{Vi v)] = ^{ x ~y)i an d by implementing secondary quan- 
tization in the so-called Fock representation. After convenient Bogoliubov transformations, 
one obtains the transition amplitudes for the vacuum state and the associated spectrum of 
the produced particles in a non-stationary background Q, Q. Many conceptual aspects 
related to possible ambiguities arising from the dependence of the particle number density 



with the reference frame and their connections with moving detectors have also attracted 
attention in the literature (see, for instance, Refs. Q, Q, 19]). 



III. EXPANDING UNIVERSE IN f(R) THEORIES 

Let us briefly recall how the evolution of the scale factor a{rf) is modified in the context 
of f{R) theories. Following Palatini's approach, where the metric and the connection are 
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independent variables, the generalized Einstein equations for a self-gravitating fluid are given 



by 



18 



2d, 



21 



28] 



f, R R^T) - i f(R) 9ia , = 8nG , 



(5) 



where f t R = df(R)/dR, R^iY) is the Ricci tensor of any independent torsionless connection 
r, R = g al3 R a f3{T) is the generalized Ricci scalar (which reduces to the ordinary Ricci scalar 
R only for some special cases) and T^ v is the energy-momentum tensor. The trace of the 
above equation reads 

f >R (R)R-2f(R)=87rGr, (6) 

where r = trT^, and for a perfect fluid it is r = p — 3p. The above identity controls solutions 
of Eq. (J5j) in the sense that, given a / function we can obtain solutions R depending on 
t. In certain special cases these solutions can be obtained explicitly. In particular, for a 
radiation or vacuum dominated universe the solution is simply R = or R = const, so that 
the solutions of the generalized Friedman equation can be explicitly obtained. 
Without loss of generality, we can choose f(R) in the following form: 



f(R)=R + g(R). 



(7) 



In the general case of a medium with equation of state in the form, p = up, one obtains a 
generalized FRW differential equation in terms of R(T). With the help of the trace identity 
as given by Eq. ([6]), we can write the following equation for the expansion of the scale factor 
in the conformal time r\: 

2 

+ a 2 G(A,R) = 0, 



+ (A-r 



(8) 



with 



G(A, R) 



9 
6 



A 
~6 



(Rg >R -g) + (A-l) 



tt2 r ' 1 

H g,R H — g,RR + 

a 4 



R!_ 

a 



2 n 2 

y,RR 



where A = (1 + 3u)/2 and H is the Hubble parameter. Note that if g(R) = then 
G(A, R) = and the general relativistic FRW differential equation for the scale factor is 



recovered 



301 ] . It should also be recalled that only in the case r = (radiation), we have 



R equal to the ordinary Ricci scalar R. As we shall see in the next section, the solutions of 
Eq. (jSJ) in this case are easily obtained for some special expressions of g(R). 
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IV. PARTICLE CREATION IN f(R) GRAVITY 



To begin with we consider the particle creation process in a radiation dominated universe, 
10 — 1/3 (or A = 1) for which R = R. In this case, Eq. (jSj) is simplified to: 



where 



— + a 2 G(l,R) = 
a 



G{l,R) = - 9 - + ±{Rg, R -g) 



(9) 



(10) 



and R = —6a" /a 3 . Inserting Eq. ([9]) into Eq. ([3]) we see that G acts exactly as an effective 
mass, and, therefore, it works like a source of particle production even in the massless case. 
Let us now assume that the g(R) function is a power law 



g(R) = f3K' 



(0 < n < 2) 



(11) 



where (3 is a dimensional constant. In this case Eq. (Q takes the form 



*" 1, 



^- + l[{3(n-2)]-—a 2 = 



or, equivalently, 



a" = —a 3 



a 



(12) 



(13) 



For simplicity, in what follows it will be assumed that (3 is positive. An interesting 
particular solution of the above equation is 



a{rj) 



VV2 



a 



V 



-OO < T] < 



(14) 



3l|. 



which has exactly the same form of a de Sitter expansion in general relativity 

Now, inserting the above result into Eqs. ([3]) and (j4j), we find that the mode function 
satisfies 



k - (l-Qam 2 ) — 
i] 2 



Xk = 



(15) 



This is a second order differential equation whose independent solutions, Xife an d X2k, are 
given by 

Xik(v) = c iVW\Jv( k \v\) , (16) 
X2M = C2^W\YMV\) , (17) 
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where J u , Y v are the Bessel functions of first and second kind, respectively, and the mass 
dependent parameter, v = a/9/4 — 12am 2 , defining the order of the functions. For all times, 
these modes must be normalized according to 

W*{ri) = X1MX2M ~ x[ k (v)X2k(v) = -2z. (18) 



It is also worth noticing that Eq. ( 1151) has some interesting features. First, it has exactly 
the same form obeyed by the mode functions in the case of a de Sitter cosmology in the 
framework of general relativity. In other words, we start with a radiation dominated universe 
in a f(R) theory and the evolution is described by the same equation of a pure de Sitter 
universe in the linear case. Particle creation in the de Sitter universe in the context of GR 



has been studied by several authors 



32 



33j, |34(. Another interesting feature of Eq. ([15!) is 



that it does not depends on the parameters (3 and n in the masssless limit. 
The Bogoliubov coefficients can be calculated and a straightforward ca 



the final expression for the number of particles created in the k mode 



32 



culation leads to 
3J: 



W-^lxMI 2 + ^lx^)l ! -i. (19) 

This expression correctly reproduces the initial vacuum condition Nk(f] — > —00) = 0. 

The concentration of created particles is readily obtained by integrating over all the 
modes: 

n= — ^ / N k ( V )d 3 k. (20) 
Now it is easy to calculate the particle creation, at least for some specific cases. 

A. Massive Case 

In order to illustrate the phenomenon of particle creation in FRW type cosmological 
models inspired by f(R) gravity, let us first consider some simple example for massive 
particles. In the case v — 1, that is, for m 2 = 5/48a, it is easy to check that the normalized 
solutions are: 

Xlk[V) ~ V k^(k\r,\)J (k\r,\) - Yoik^DMklvl) ' " ' 



XM = iJl , Yim) , (22) 

V k^YMv\)Mk\v\)-Y,{k\v\)Uk\v\) 
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FIG. 1: Massive Case. Spectrum of the created particles as function of the conformal time r/ for 
two different values of the wavenumber k, from left to right, k = 1 and k = 2 respectively. The 
spikes represent the points where the frequency vanishes and iV& diverges. 

and the corresponding spectrum of created particles is readily obtained by inserting the 
above solutions into (TI9")) . 

In Fig. 1, we show the spectrum for two different wavenumbers. The spikes represent the 
points where the frequency vanishes and diverges. 

Another interesting case is m 2 = I /6a which is conformally related to flat space-time. 
The frequency vanishes for this choice of mass so that the modes are plane waves at all 
times, and, therefore, there is no particle production. The absence of matter creation in this 
f{R) massive case is similar to the classical result derived by Parker in general relativity 
for the massless case [ijj]. Implicitly, it also suggests the possibility of massless particles 
production during the radiation dominated phase in f(R) models. In order to show that, 
let us now calculate the spectrum of the created massless particles. 
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FIG. 2: Massless Case. Spectrum of created particles as a function of the conformal time rj for 
different values of the wavenumber k, from left to the right, k = 1, k = 2 and k = 3 respectively. 
The spikes represents the points where the frequency vanishes and diverges. 

B. Massless Case 

The main result of this paper is related to the creation of massless particles (m = 0). In 
this case, the normalized mode functions (I16p and (JT7J) takes the following forms 

xM =7t( 1+ i) e ' tv ' (23) 
^ = 7e - £) (24) 

Now, inserting the above expressions into Eq. (fT9l) we obtain the spectrum of created 
particles for each mode k: 

Nk(v) = ~2 - i'^V " 2 ' 1/2 (foj + jfc^s) ~ 4 |A;V-2| 1 /2 ( kr > ~ ^ + jfc^a) ' 

(25) 

In Fig. 2 we display the graphic of this function for some values of k. Apart from the 
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spikes which occur due to the vanishing of the frequency, we also see that the number of 
created particles grows in the limit i] ^ 0. 

V. CONCLUDING REMARKS 

Due to the discovery of the accelerating Universe at low redshifts, an increasing attention 
has been paid to f(R) gravity as a possibility to reduce the so-called cosmological dark 
sector. In this kind of cosmology the Universe can be accelerating driven only by cold dark 
matter. 

In this work, assuming that the deviations from general relativity are described by a 
power law, we have addressed the problem of particle creation in f(R) homogeneous and 
isotropic cosmologies dominated by a fluid obeying a cj-type equation of state. The mode 
equation in such background was derived and their general solution explicitly obtained. 
As an illustration, we have calculated the spectrum of created particles for a radiation 
dominated universe considering the cases of massive and massless particles. 

Unlike the Parker results which are based on the standard general relativity, we have also 
explicitly shown that massless particles can be produced during the radiation dominated 
phase in the framework of f{R) cosmologies. This result is not valid only for the case of a 
scalar field as discussed here. Actually, due to the new contributions to the FRW differential 
equation, one may show that the positive- and negative-frequency parts of the fields become 
mixed in an expanding f(R) cosmology thereby leading to the creation of massless particles 
of nonzero spin even at early times. 

Acknowledgments 

JASL is partially supported by CNPq and FAPESP No. 04/13668-0 and CHGB is sup- 
ported by CNPq No 150429/2009-6 (Brazilian Research Agency). 



[1] R. Brout et. al. Ann. Phys. (N.Y.), 115, 78 (1978); D. Atkatz and H. Pagels, Phys. Rev. D 25, 
2065 (1982); J. B. Hartle, in The Very early Universe, edited by G. Gibbons, S. W. Hawking, 
and S. Siklos (Cambridge UP, Cambridge, 1983); B. L. Hu and Pavon, Phys. Lett. B 180, 329 
(1986). 

[2] Ya. B. Zeldovich, JETP Lett. 12, 307 (1970). 

11 



[3] G. L. Murphy, Phys. Rev. D48, 4231 (1973); Z. Klimek, Acta Cosmologica 3, 49 (1975). 
[4] B. L. Hu, Phys. Lett. A90, 375 (1982). Adv. Astrop. 1, 23 (1983). 

[5] J. D. Barrow, Phys. Lett B180, 335 (1986); M. Morikawa and M. Sasaki, Phys. Lett. B 165, 

59 (1985); T. Padmanabhan and S. M. Chitre, Phys. Lett. A 120 433 (1987); J. A. S. Lima, 

R. Portugal and I. Waga, Phys. Rev. D 37, 2755 (1988). 
[6] I. Prigogine, J. Geheniau, E. Gunzig, P. Nardone, Gen. Rel. Grav. 21, 767 (1989). 
[7] M. O. Calvao, J. A. S. Lima, I. Waga, Phys. Lett. A162, 223 (1992). See also, J. A. S. 

Lima, M. O. Calvao, I. Waga, "Cosmology, Thermodynamics and Matter Creation", Frontier 

Physics, Essays in Honor of J ay me Tiomno, Edited by , (World Scientific, Singapore, 1991). 
[8] J. A. S. Lima, A. S. M. Germano, Phys. Lett. A170, 373 (1992); W. Zimdahl and D. Pavon, 

Mon. Not. R. Astr. Soc. 266, 872 (1994); W. Zimdahl and D. Pavon, Gen. Rel. Grav. 26, 

1259 (1994); L. R. W. Abramo and J. A. S. Lima, Class. Quant. Grav. 13, 2953 (1996). 
[9] J. Gariel, G. Le Denmat, Phys. Lett. A200, 11 (1995); J. A. S. Lima, A. S. M. Germano and 

L. R. W. Abramo, Phys. Rev. D53, 4287 (1996). 
[10] J. A. S. Lima and J. S. Alcaniz, Astron. Astrophys. 348, 1 (1999), [astro-ph/9902337] ; Zimd- 

hal, D. J. Schwarz, A. B. Balakin and D. Pavon, Phys. Rev. D 64, 063501 (2001). 
[11] J. A. S. Lima, F. E. Silva and R. C. Santos, Class. Quant. Grav. 25, 205006 (2008), 

larXiv:0807.3379l [astro-ph] ; G. Steigman, R. C. Santos and J. A. S. Lima, JCAP 06033 (2009), 

larXiv:0812.39"T2l [astro-ph]. 
[12] L. Parker, Phys. Rev. Lett. 21, 562 (1968). 

[13] L. Parker, Phys. Rev. 183, 1057 (1969); Phys. Rev. D 3, 346 (1971); Phys. Rev. Lett. 28, 

705 (1972); Phys. Rev. D 7, 976 (1973). 
[14] N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space (Cambridge University 

Press, Cambridge, 1982); S. A. Fulling, Aspects of Quantum Field Theory in Curved Spacetime 

(Cambridge University Press, Cambridge, 1989); A. A. Grib, S. G. Mamayev and V. M. 

Mostepanenko, Vaccum Quantum effects in Strong Fields (Friedmann Laboratory Publishing, 

St. Petesburg, 1994). 

[15] V. F. Mukhanov and S. Winitzki, Introduction to Quantum Effects in Gravity, (Cambridge 

University Press, Cambridge, 2007). 
[16] Ya. B ZePdovich, Pisma Zh. Eksp. Teor. Fiz. 12, 443 (1970), (English transl. JETP 12, 307 

(1970)); A. A. Grib, S. G. Mamayev and V. M. Mostepanenko, Gen. Rel. Grav., 7, 535 (1975); 

12 



A. A. Grib and Yu. V. Pavlov, [gr-qc/0505140] ; Grav. Cosmol. 11, 119 (2005); Grav. Cosmol. 
12, 159 (2006). 

[17] L. P. Grishchuk, Class. Quant. Grav. 10, 2449 (1993); M. R. G. Maia, Phys. Rev. D 48, 647 

(1993); M. R. G. Maia and J. D. Barrow, Phys. Rev. D 50, 6262 (1994); M. R. G. Maia and 

J. A. S. Lima, Phys. Rev. D 54, 6111 (1996). 
[18] A. A. Grib and S. G. Mamayev, Yad. Fiz. 10, 1276 (1969) [English transl.: Sov. J. Nucl. Phys. 

10, 722 (1970)]. See also Yu. V. Pavlov, Theor. Math. Phys. 126, 92 (2001), [gr-qc/0012082| . 
[19] M. I. Shirokov, Sov. J. Nucl. Phys. (USA), 7, 411 (1968); A. A. Grib and S. G. Mamayev, 

Sov. J. Nucl. Phys. (USA), 14, 450 (1972); J. R. Letaw and J. D. Pfautsch, J. Math. Phys. 

23, 425 (1982); W. Junker and E. Schrobe, Ann. Inst. Henry Poincare 3, 1113 (2002); L. C. 

Crispino, A. Higuchi and G. E. A. Matsas, Rev. Mod. Phys. 80, 787 (2009). 
[20] D. N. Vollick, Phys. Rev. D 68, 063510 (2003); M. Amarzguioui, O. Elgaray, D. F. Mota and 

T. Multamaki, Astron. Astrophys. 454, 707, (2006); S. Nojiri and S. D. Odintsov, Phys. Lett. 

B 659, 821, (2008). 

[21] G. Allemandi, A. Borowiec, M. Francaviglia and S. D. Odintsov, Phys. Rev. D 72, 063505 
(2005); L. Amendola, D. Polarski and S. Tsujikawa, Phys. Rev. Lett. 98, 131302 (2007); J. 
Santos, J. S. Alcaniz, F. C. Carvalho and N. Pires, Phys. Lett. B 669, 14 (2008); J. Santos 
and M. J. Reboucas, Phys. Rev. D 80, 063009 (2009). For a review see, T. P. Sotiriou and V. 
Faraoni. [arXiv:0805.172"6V 2. [gr-qc], Rev. Mod. Phys. (2009). To appear. 

[22] M. Gasperine and G. Veneziano, Phys. Lett. B 277, 256 (1992). 

[23] G. A. Vilkovisky, Class. Quant. Grav. 9, 895 (1992). 

[24] A. A. Starobinsky, Phys. Lett. B 91, 99 (1980). 

[25] A. D. Dolgov and M. Kawasaki, Phys. Lett. B 573, 1 (2003); T. Chiba, Phys. Lett. B 575, 
1 (2003). 

[26] J. A. R. Cembranos, Phys. Rev. D 73, 064029 (2006); T. P. Sotiriou, Class. Quant. Grav. 23, 
1253 (2006). 

[27] T. P. Sotiriou, Phys. Rev. D 73 063515 (2006). 

[28] X. Meng and P. Wang, Class. Quant. Grav. 21, 2029 (2004). 

[29] V. Miranda, S. E. Joras, I. Waga and M. Quartin, Phys. Rev. Lett. 102, 221101 (2009). 
[30] M. J. D. Assad and J. A. S. Lima, Gen. Rel. Grav. 20, 527 (1988); J. A. S. Lima, Am. J. 
Phys. 69, 1245 (2001). 

13 



[31] V. Mukhanov, Physical Foundations of Cosmology, Cambridge University Press, (2005). 
[32] M. Bordag, J. Lindig and V. M. Mostepanenko, Class. Quantum Grav. 15 581 (1998). 
[33] E. Mottola, Phys. Rev. D 31, 754 (1985); C. Molina-Pans, Int. J. Th. Phys. 38, 1273 (1999). 
[34] M. Mijic, Phys. Rev. D 57, 2138 (1998); S. Biswas and I. Chowdhury, Int. J. Mod. Phys. D 
15, 937 (2006). 



14 



